Metal-insulator transition in disordered 2DEG 
including temperature effects. 
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We calculate self-consistently the mutual dependence of electron correlations and electron-defect 
scattering for a two dimensional electron gas at finite temperature. We employ an STLS approach 
to calculate the electron correlations while the electron scattering rate off Coulombic impurities 
and surface roughness is calculated using self-consistent current-relaxation theory. The methods 
are combined and self-consistently solved. We discuss a metal-insulator transition for a range of 
disorder levels and electron densities. Our results are in good agreement with recent experimental 
observations. 
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In recent experiments on two-dimensional elec- 

tronic systems in zero magnetic field a well denned metal- 
insulator transition has been observed. The transition 
contradicts the prediction pj that two dimensional elec- 
tron systems are always localized in the presence of any 
disorder. The mobility is large at the transition point 
(fi > lm 2 /Vs). For such a high mobility the scattering 
off defects will be weak and electron correlations could 
make a significant difference. 

Both disorder and correlations by themselves can lead 
to at least two different types of localization. In the pres- 
ence of disorder Abrahams et al Q showed that non- 
interacting electrons in two dimensions cannot sustain 
static conductivity no matter how small the level of dis- 
order is. On the other hand if there are strong correla- 
tions between the charge carriers this can cause a differ- 
ent type of localization in disorder-free systems which 
is associated with Wigner crystallization 0. In real 
systems there are both electron-electron and electron- 
impurity interactions. With weak disorder and in the 
low electron density limit the electron-electron interac- 
tions should dominate leading to Wigner crystallization, 
while in the high electron density limit the electron cor- 
relations are very weak and the localization should be of 
the Anderson type. Between these two extremes disorder 
and correlations compete with each other to decide the 
nature of the localization. 

In this paper we examine the interdependence of corre- 
lations and defect scattering at finite temperatures. We 
use the self-consistent formalism of Singwi, Tosi, Land 
and Sjolander (STLS) || to treat electron-electron cor- 
relations, while for the electron-defect scattering we use a 
memory function approach [|7| to calculate the decay time 
of the density fluctuations from scatterings off defects. 
Previous calculations of the influence of disorder have ei- 
ther used the first Born approximation or else they have 
introduced an adjustable parameter 7 for the electron 
scattering rate from the disorder. In these approaches 
the scattering rate is not affected by correlations while 
the two effects are in fact interdependent and should be 



self-consistently linked. 

We consider a random distribution of Coulombic im- 
purities with charge e and density Uj. The charge car- 
rier (electron or hole) density n c is chosen low enough 
that the carriers occupy only the lowest energy sub- 
band. In the STLS formalism the bare Coulomb potential 
V{q) = 2ire 2 /qe acting between the carriers is replaced by 
an effective interaction V e g(q) = V(q)[l — G(q)], where 
G(q) is a local field factor. V e g(q) is generally weaker 
than V(q) because of the correlations between the carri- 
ers. The response function becomes 
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We include the effect of disorder by replacing the free 
particle dynamical susceptibility x' '' (?> w ) m Eq. |l| by 
the response function for non-interacting charge carriers 
scattering from the disorder B , 
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X (?) <*> + 17) is the susceptibility for non-interacting 
carriers scattering off the disorder. We use the mem- 
ory function formalism j9) to calculate the scattering rate 
7. This is related to the mobility /i by the Drude ex- 
pression fj, = e/m*7. At zero temperature in the dif- 



fusive regime x (<7j w ) is given by lim a 
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(2m*)/(TTk F n 2 )(Vq 2 )/(Vq 2 + iuj), where V = v 2 F /-f is the 
diffusion constant. In the limit when 7 goes to infin- 
ity the system becomes non-diffusive. This represents a 
localized phase Q . 

In the memory function formalism 7 is expressed in 
terms of the non-linear equation H 

i7 = -r^E? 2 Wimp(q)\ 2 ) + (|^ SUI f(g)| 2 )] 

) 711 * 11 _ ' L J 



m 

x (c %) 



$(q,i 7 ) 



(3) 



1 



<f>{q,i-y) 



0o (gihO 

l + i70o(g,i7)/x (o) (9) 



(4) 



< &(g, i 7 ) is the relaxation spectrum for non-interacting 
carriers that scatter off the disorder with scattering 
rate 7. It is expressed in terms of 0o(<Li7) = 
(I/17) [x^ (?) 17) — X (q)] which is the relaxation func- 
tion for non-interacting carriers with lifetime 7 . 
The static response function x{q) — x^(q)/{l + 
— G{q)}x^°Hq)} includes the correlations between 
the carriers for the disordered system. Ui mp (q) — 
[(27re 2 )/(eg)] exp(— qd)Fi(q) is the impurity potential for 
the monovalent Coulomb impurities which are in a layer 
separated from the electron or hole plane by a distance d. 
We use for the impurity form factor Fi(q) Eq. 4.28 in Ref. 
jn|. For electrons in the Si MOSFETs we also include 
interface surface roughness scattering. This is the term 
W sm{ (q) = % /7rAAr(g)exp(-(qA) 2 /8). Values for the pa- 
rameters A = 0.37nm and A = 2.0nm are taken from Si 
MOSFET data @. For T{q) we use the expression in 
Ref. . For GaAs surface roughness scattering is much 
smaller and we set W sur f(q) = 0. 

In the STLS formalism the density-density correla- 
tion function (6n(r,t)6h(r' ,£)) is approximated by the 
non-linear product, 6n(r,t)xg(r — r')x5n(r' ,t). The 
5n(r,t) are expectation values and g(r) is the pair- 
correlation function giving the probability of finding 
electrons a distance r apart. Using g(r) = 1 + 
n^ 1 J d 2 q exp(iq.r)[S'(g) — 1], this gives us a relation be- 
tween the static structure factor S(q) and the local field 
factor G(q), 
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1 f d 2 k (q-k)F(k) 
(2tt) 2 q 2 V(q) 



[5(|q-k|)-l]. (5) 



For a given G(q) we can determine 7 from Eq. |^. Then 
S(q) = (n c 7r) _1 J* °° did lmx{q,uj) is calculated from the 
fluctuation-dissipation theorem and 



x(q,u) 



X {s) (q,u) 
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Equation || gives a new local field factor which can be 
used in Eq. H. The process is repeated iteratively until 
there is overall self-consistency with both 7 and G(q). 

Our variables are the carrier density n c , the impurity 
density n,, and the distance separating the impurities 
from the electron or hole plane d. The Metal Insulator 
Transition is observed at carrier densities that are rel- 
atively low by conventional semiconductor standards so 
the Fermi temperature TV = -Ef/^b can be of the order 
of a few degrees K. We have solved the equations at fi- 
nite temperatures and determined the dependence of 7 
on temperature up to T < Tp. 




FIG. 1. Local-field-factor G(q) for different impurity den- 
sities. The impurities are separated from the carrier plane by 
d=5ag. Carrier density is n c = 35 x 10 10 cm -2 . 

Figure |l| shows the dependence of the local field fac- 
tor G(q) on the impurity concentration at zero temper- 
ature. Increasing the disorder enhances G{q). This is 
caused by the decrease in x as the scattering rate 7 
gets bigger (Eq. ||). Enhancing G(q) weakens the effec- 
tive interaction between the carriers, and hence weakens 
the screening of the carrier-impurity potential. The net 
result of enhancing G(q) is thus to strengthen the effect 
of the disorder potential. This in turn further increases 
7. At a critical level of disorder this non-linear feedback 
causes 7 to increase rapidly and leads to localization of 
the carriers. 
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FIG. 2. Scattering rate 7 for Si as a function of electron 
density for impurity density 0.5 xlO 11 cm -2 . Separation of 
the impurity layer is d=a,Q. Curves are for temperatures T 
(labels are in units of the Fermi temperature Tf for density 



We see this non-linear behavior in Fig. At a cer- 
tain critical carrier density the scattering rate 7 starts 
to increase rapidly. The impurity density here is n% = 



2 



0.5 x lC^cm -2 , with the impurities separated from the 
carrier plane by distance d = o,g (aj is the effective Bohr 
radius) . The non- linear increase in 7 is due (i) to the en- 
hancement of the local field G(q) which strengthens the 
disorder potential, and (ii) to the rapid increase in the 
relaxation spectrum < J > (g, vy) with 7. Figure || also shows 
the dependence of 7 on temperature. The labels on the 
curves give the temperature in units of Tp for density n c . 
7 diminishes with increasing temperature. This reflects 
the weakening of the correlations at finite T. 
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FIG. 3. Scattering rate 7 for Si at T = as a function 
of electron density. Curves are labelled for different impu- 
rity densities m (in units of 10 11 cm -2 ). Impurity separation 

In Figure || we show the increase of 7 with impurity 
density at T — 0. The curve labels give rij in units of 
10 11 cm -2 . The impurity separation is d = 5<2q. As ex- 
pected, increasing the impurity concentration has the ef- 
fect of increasing the scattering rate. On the other hand, 
increasing the carrier density has the opposite effect on 7 
because the correlations are reduced. However the non- 
linear effects from n c are much stronger than the non- 
linear effects from n.^. 

There exists no first order localization-delocalization 
transition in two-dimensions. In the metallic phase the 
mean- free-path is large £kp ^> 1 while in the localized 
phase £kF should be < 1. We take the localization 
boundary to be the point where Ikp = 1. Table I gives 
for different temperatures the critical electron density at 
which ikp passes through unity. When the system is at 
finite temperature the correlations are reduced and one 
needs to go to slightly lower carrier densities before lo- 
calization can be achieved. 

Table I shows the temperature dependence is quite 
small and so we give Metal-Insulator phase diagrams 
for T = 0. The dependence of the hole density at the 
transition on the impurity layer separation d for GaAs 
is shown in Fig. m The impurity density is fixed at 
rii = 0.5xl0 11 cm . We compare our results with the 



observation of the transition in Ref. j| in a GaAs sam- 
ple. Our curve is consistent with this measurement. The 
phase boundary is sensitive to d because of the exponen- 
tial factor in Ui mp (q). This suppresses carrier-impurity 
scattering for short wavelengths q^d. For the residual 
long wavelength scattering the correlations are always 
weak. Thus when d is large the correlations play a rela- 
tively less important role in the localization. 
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FIG. 4. Phase diagram for holes in GaAs. Critical hole 
density at which 7 = 1 as a function of impurity layer 
separation d. Impurity density is n; = 22xl0 11 cm -2 . The 
experimental data point for GaAs is from Ref. M. 
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FIG. 5. Phase diagram for electrons in Si. Critical carrier 
density at which 7 = 1 as a function of impurity density for 
separations d=0 and 5oq. The experimental data points are 
taken from Ref. Q|. 

The phase diagram in Figure ^| plots the critical elec- 
tron density in Si at the transition as a function of im- 
purity density n^. The impurity layer separations are 
d = and 5oq. We also show experimental data points 
for the position of the Metal-Insulator transition in Si 
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0). Again our predicted phase boundary agrees with 
these observations. On our curve the electron-impurity 
scattering becomes weaker as the critical electron den- 
sity decreases. Thus at the transition point the electron 
mobility will increase as the electron density decreases. 
This is consistent with recent observations g. 

We conclude that correlations and impurity scattering 
do mutually affect the localization transition. Finite tem- 
peratures tend to suppress correlations and this slightly 
reduces the critical carrier density for localization to oc- 
cur by an amount which becomes significant for temper- 
atures of the order Tp. 
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